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The PHOBOS experiment at RHIC has measured the multiplicity of primary charged 
particles as a function of centrality and pseudorapidity in Au+Au collisions at ^s^n = 
19.6, 130 and 200 GeV. Two kinds of universal behavior are observed in charged particle 
production in heavy ion collisions. The first is that forward particle production, over a 
range of energies, follows a universal limiting curve with a non-trivial centrality depen- 
dence. The second arises from comparisons with pp/pp and e + e~ data. (N ch ) / (N part /2) in 
nuclear collisions at high energy scales with y^i in a similar way as N c h in e + e~ collisions 
and has a very weak centrality dependence. This feature may be related to a reduction in 
the leading particle effect due to the multiple collisions suffered per participant in heavy 
ion collisions. 

1. Introduction 

The PHOBOS experiment has measured dN c h/dr) and the average multiplicity of charged 
particles (N c h) produced in heavy ion collisions for center of mass energies in the nucleon- 
nucleon center of mass system, ^snn, of 19.6, 130 and 200 GeV. The data is also binned 
as a function of event centrality (impact parameter) characterized by the number of par- 
ticipating nucleons, N part , allowing comparisons to elementary systems, like pp/pp and 
e + e~ — > hadrons. 
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Figure 1. dN/dr/ / (N part /2) for peripheral Figure 2. Particle density at midrapidity 

and central events at three RHIC energies. t A . A /_ j + - 

to lor A+A, pp/pp and e^e . 

The PHOBOS multiplicity detector consists of several arrays of silicon detectors which 
cover nearly the full solid angle for collision events. The "Octagon" detector surrounds the 
interaction region with a cylindrical geometry below \z\ < 50 cm, covering \r)\ < 3.2. Two 
sets of 3 "Ring" detectors are placed far forward and backward of the interaction point 
and surround the beam pipe, covering 3 < \rj\ < 5.4. The methods used for measuring 
the multiplicity of charged particles as well as for extracting (N part ) has been described 
in more detail in Ref . [ [1|] . 



2. Limiting Behavior in Pseudorapidity Distributions 

Fig. [I] shows dN c h/dr}' / {N par t /2) (V = V — 2/beam) measured at three different RHIC 
energies for peripheral ({N part ) ~ 100) and central events {{N part ) ~ 355), in the left and 
right panels, respectively. These show a clear "limiting behavior" in the fragmentation 
region. That is, the distributions are independent of beam energy over a substantial range 
in t]'. As the beam energy increases, dN/drj' follows the universal trend until it reaches 
85-90% of its maximum value at midrapidity, at which point it stops following the trend. 
Similar behavior has been observed in elementary collisions as well, both in pp collisions 
[ 0] and in e + e~ annihilation to hadrons [ |J. 

The limiting curve constrains the energy dependence of the charged particle multiplicity. 
It also varies with centrality in such a way that the increases seen at low (which is 
midrapidity in rj) as N part increases, are accompanied by decreases near rj' ~ (forward 
rapidities), as seen in Fig. |l|. It is not clear why this behavior occurs, e.g. from energy 
conservation or some kind of long-range correlation. 



3. Comparison with Elementary Systems 

Comparisons of the plateau height, dN c h/ dr]^^/ {Npart/2), in heavy ion collisions with 
pp data have been made previously [ |4[]. However, we can also include data from e + e~, 
presented as dN/dyx, the rapidity density along the event thrust axis, calculated assuming 
the pion mass. JETSET calculations suggest that the difference between dN/dyx and 
dN/dr] is less than ±10% for \y\ < 4. In the comparison at midrapidity, shown in Fig. 
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Figure 3. Top: dN c h/dr] for central 
Au+Au and collisions compared with 
dN/dy T for e + e~ data, all at ^/s^n = 200 
GeV. Bottom: Au+Au and pp data di- 
vided by a Woods-Saxon fit to the Au+Au 
data. 
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Figure 4. Comparison of (AT c/l )/ (N part /2) 
for A+A, pp/pp, and e + e~ data compared 
with a fit to the e + e~ data. 



^ we find two interesting features. First, the energy dependence of all of the systems is 
approximately logarithmic, at least below 200 GeV. Secondly, while it has been noticed 
that heavy ion data is 40-50% above pp/pp data, the e + e" data has the same trend and 
a similar magnitude (within 10%) as AA over a large range in ^fs = 14 — 183 GeV [ 0. 
This correspondence holds over the bulk of the angular distribution, as shown in Fig. |3|, 
where central Au+Au (divided by (N part /2)), pp [0 and e + e _ [|J data are compared, all 
at y/s = 200 GeV. In the lower panel, we observe that the shapes of the Au+Au and pp 
distributions are also very similar over a large r] range, but the integrals differ by ~ 40%. 

In Fig. £|, we compare (N c h) / (N part /2) in heavy ion collisions [ § to e + e _ and pp/pp 
data over a large range in y/s [0. It is observed that (N c h) / (N part /2) lies below pp at low 
energies, passes through the pp data around yfs ~ 10 GeV, and then gradually joins with 
the e + e~ trend above CERN SPS energies. These comparisons can be seen more clearly 
by dividing all of the data by a fit to the e + e~ data [ The pp/pp data follows the same 
trend as e + e~, but it can be shown that it matches very well if the "effective energy" 
V^eff = is used, which accounts for the leading particle effect seen in pp collisions [ 

Ref. [ finds that bulk particle production in pp and e + e~ data does not depend in 
detail on the collision system but rather the energy available for particle production. In 
this scenario, the Au+Au data suggests a substantially reduced leading particle effect in 
central collisions of heavy nuclei at high energy. 

The alleviation of the leading particle effect might not be so surprising in nuclear 
collisions. Each participating nucleon is struck V > 3 times on average as it passes 
through the oncoming gold nucleus for N part > 65. One could speculate that the multiple 
collisions transfer much more of the initial longitudinal energy into particle production. 
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This naturally leads to the scaling of total particle production in heavy ion collisions with 
Npart, as seen i n Fig- B reminiscent of the "wounded nucleon model" [ [10] but with the 
scaling factor determined by e + e~ rather than pp. 



4. Conclusions 



In conclusion, PHOBOS has observed 
two kinds of universal behavior. The first 
is an energy-independent, but centrality- 
dependent, universal limiting distribution 
of charged particle production away from 
midrapidity. The second is that the to- 
tal charged particle multiplicity per par- 
ticipant pair in heavy ion collisions above 
CERN SPS energies scales with yfs in a 
similar way as e + e~ collisions. These two 
kinds of universality strongly constrain the 
energy dependence of the total charged par- 
ticle multiplicity and angular distributions, 
and may offer a new perspective on parti- 
cle production in heavy ion collisions. 
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Figure 5: The total number of charged par- 
ticles per participant pair shown as a func- 
tion of N part for y/SNN = 19.6, 130, and 200 
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